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Abstract 

Quantum phase transitions occur when the ground state of a quantum system undergoes a qual- 
itative change when an external control parameter reaches a critical value. Here, we demonstrate a 
technique for studying quantum systems undergoing a phase transition by coupling the system to a 
probe qubit. It uses directly the increased sensibility of the quantum system to perturbations when 
it is close to a critical point. Using an NMR quantum simulator, we demonstrate this measurement 
technique for two different types of quantum phase transitions in an Ising spin chain. 
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Introduction. — Phase transitions describe sudden changes in the properties of a physi- 
cal system when an external control parameter changes through some critical value. If the 
system under consideration is a quantum mechanical system in its ground state, i.e. at 
zero temperature, and the phase transition occurs as a function of a non-thermal control 
parameter, we speak of quantum phase transitions (QPTs) Examples include the tran- 



sitions in superconductors 



and fractional quantum Hall systems l3|] . Related phenomena 



a. 



have also been experimentally observed in heavy Fermion systems [^, common metals 
and in Bose- Einstein condensates 6|. QPTs occur as a result of competing interactions 
and the different phases often show different types of correlations between the constituents, 
with correlation lengths that can become arbitrarily large. When specific quantum effects 
of phase transitions are of interest, it is therefore natural to compare entanglement in the 
different phases j?, 8]. 

Experimental observations of QPTs are relatively straightforward when they are accom- 
panied by a change of a suitable order parameter, such as the conductivity or susceptibility 
in superconductors or the total magnetization in some spin chains j3, S]- However, such 
global measurements cannot provide all the details and they are not suitable for closer in- 
vestigations of the systems in the interesting area close to the critical points. Moreover, 
not all order parameters can be measured by global measurements. A complete analysis of 
the system is provided by quantum state tomography j^, [l^ , but this approach scales very 
poorly with the size of the system. 

As a possible alternative for closer investigations of quantum systems in the vicinity of 
critical points, it was suggested to compare the evolution of systems at slightly different 
values of the control parameter. This approach may be considered as a visualization of 
" quantum fluctuations" . Different possibilities exist for compa ring these evolutions, some of 
which have been called Loschmidt echo (LE) or fidelity decay In the vicinity of critical 
points, the systems are expected to be much more susceptible to external perturbations than 
in the center of a phase 12]. Such a comparison is possible by coupling the system under 
study to a second quantum system, consisting in the simplest case of a single qubit. The two 
states of the probe qubit can then be used to probe the system under two different values 
of the control parameter. The signal obtained in this case corresponds to the overlap of two 
states evolving under slightly different control parameters. 

In this Letter we implement this protocol in a nuclear magnetic resonance (NMR) quan- 
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turn information processor. The system undergoing the QPT corresponds to an Ising-type 
spin chain and the control parameter to a longitudinal magnetic field. In a purely longitu- 
dinal field, the ground state is degenerate at the critical points. This degeneracy is lifted if 
the magnetic field contains a transverse component. For the longitudinal as well as for the 
transverse case, we measure the QPT by coupling the spin chain to a probe qubit. 

Level-crossing. — We first consider a QPT in the Ising model in a minimal system con- 
sisting of two spins 1/2. Its Hamiltonian is 

H^ = alal + B,ial + a',), (1) 

where the a* are Pauli operators and Bz is a magnetic field. The units have been chosen 
such that the coupling constant between the two qubits is 1. For the purpose of this paper, 
it is sufficient to consider the triplet manyfold. Within this subsystem, the ground state 
depends on the field strength: 

' |00) (5, < -1) 
1^,(5.)) = 10+) (-1 <Bz<l) (2) 
. Ill) (Bz > 1), 

where = (|01) + \10))/y/2. Figure [D shows the energy levels of the system and the 
concurrence of the ground state. Obviously B^ = ±1 are the critical points. The low-field 
phase is maximally entangled (C = 1), while the high-field phases correspond to product 
states (C = 0). 

In this system, the QPTs occur at points where the ground state is degenerate. Close 
to this critical point, it is therefore very susceptible to small perturbations. If we couple it 
to a probe qubit (which we label 0) via the interaction ecr^{al + al), the total (three-qubit) 
system can be decomposed into two subsystems, in which qubits 1 and 2 "see" an effective 
ield Bz ±5. If these two fields fall on different sides of the critical point, the "state overlap" 
3] L = \ {il)g{Bzfl)\il)g{Bz^i))\^ vanishes, otherwise it is unity, as shown by the thin line in 
Figure [H Here, Bz^ = Bz + e specifies the effective field for the subsystem coupled to |0)o 
and correspondingly for the other subsystem. In the extreme case where the two states 
of the probe qubit are orthogonal (L = 0), the probe qubit has "measured" the quantum 
system jl4|. 

To measure L, we first initialize the system and probe qubits into the ground state 
1 000). From there a Walsh-Hadamard transform places the probe qubit into the symmetric 



superposition state |+) = (|0) + |l))/v^. We then use the interaction between the probe and 
the system to apply a conditional evolution to the two system qubits 1 and 2: if the probe 
qubit is in state |0), the system evolves from |00) \tpg{Bzfi)), and if qubit is in state 
|1), the system evolves from |00) \tljg{Bz^i)) . The network representation of this process 



is shown in Figure [2]^a) 15|. Pq and Pi denote the conditional evolutions. The output of 
the network is 

|vl/) = [|0)|^,(S,,o)) + |l)|V^,(5.,i))]/v^. (3) 

Taking the trace over the (12)-system, one obtains the reduced density matrix p(°) of the 
probe qubit. The off-diagonal elements are = P2i''^ = {'^9{Bz,o)\4'g{Bz,i)) ■ Hence the 
overlap L can be obtained by measuring L = A\{a^)\'^, the transverse magnetization of the 
probe qubit, which can be observed as a free induction decay. 

For the experimental implementation, we chose the nuclear spins of ^^C, ^H, and ^^F 
of Diethyl-fluoromalonate as qubits, shown in Figure [2]^b). The scalar coupling constants 
are J12 = 47.6 Hz, Jiq = 161.3 Hz and J20 = —192.2 Hz. The sample consisted of a 2.3:1 
mixture of unlabeled Diethyl fluoromalonate and d6-acetone. Molecules with a ^^C nucleus, 
which we used as the quantum register, were therefore present at a concentration of about 
0.7%. 

The effective pure state |000) was prepared by spatial averaging We implemented 

the quantum network of Figure [2](a) for five cases corresponding to Bz = —1.5, —1, 0, 1, 
1.5, respectively. As an example. Figure [2](c) shows the pulse sequence when Bz = —1. 
Figure [3] shows the experimental results. The spectra on top were measured for the values 
of Bz = -1.5, 0, and +1, and the asterisks indicate the integrated signal amplitudes. Clearly, 
the integrated signal essentially vanishes at the quantum critical point, while it remains 
close to the maximum inside the three phases, in excellent agreement with the theoretical 
expectation. 

Avoided level-crossing. — If we add a transverse field to the system, the QPTs are no 
longer singular points, but they acquire a finite width. The modified Hamiltonian is 

= alal + B,{al + al) + B^ia] + a]). (4) 

As long a.s Bx ^ 1, the ground states are very close to those of Eq. ([2]), except in the 
vicinity of the critical points, where the transverse field mixes them, thus avoiding the level 
crossing. In this region, it is sufficient to consider the two lowest energy states. They form 



a two-level system that can be described by the effective Hamiltonian 

Heff - -\B,\I + (1 - \B,\)a, + V2B^a^, 

where / denotes the unit operator. Within this approximation, the ground state is 
\ijg{B^, B,)) = \U) cos(v7/2) + 10+) sin(v?/2), where tan<^ = ^25^/(1 - \B^\). When B^ > 0, 
\ll) = |11); when B, < 0, |//) = |00). 

The coupling of this system to a probe qubit provides us with a natural way of measuring 
the phase transition. As before, we use an Ising-type interaction, which results in the total 
Hamiltonian 

H^H^ + ea'^Aal + a!). (5) 

To measure the QPT, we first initialize the system into the ground state \ipg{Bx, B^)) in 
a given longitudinal field B^. As we turn on the coupling to the probe qubit initially in 1+), 
the combined system splits into two subsystems, corresponding to the two eigenstates of the 
probe qubit. In the two subsystems, the effective longitudinal field acting on qubits 1 and 2 
is 1 — \Bz\:Le, i.e. the coupling shifts the two subsystems in opposite directions along the 
Bz-axis. The eigenstates of the two subsystems arc therefore also different. In terms of the 
mixing angle ^9, the sensitivity of the basis states to the variation of the longitudinal field 
can be quantified as 

d^/d\B,\ = V2B,/[2Bl + (1 - \B,\f]. (6) 

Apparently, this is a resonant effect: The sensitivity reaches a maximum at the QPT and 
falls off with the distance from the critical points like a Lorentzian. The full width of this 
"resonance line" is equal to the splitting 2\f2Bx of the two lowest energy levels at the critical 
point. 

To measure this behavior, we initialize the probe qubit into the 1+) state. As we turn 
on the coupling to the system, each subsystem is no longer in an eigenstate, but starts to 
evolve in its new basis. The initial state \-\-)\'^g{Bx^Bz)) evolves as 

|*(r)) = [|0)|*o(r)) + |l)|*i(r))]/V2. (7) 

Here \^q{t)) describes the two system qubits coupled to the state |0) of the probe qubit, 
evolving under the Hamiltonian ifo = -f^T + ^("^l + ^V)- Similarly, \^\{t)) describes the 
probe qubit in state |1) and evolves under i^i = — e{a\ + cr^). We use this differential 
evolution for measuring the QPT via the overlap L — | (^o('?') I ('?')) P- 
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Figure H] shows the quantum circuit and pulse sequence that implement this mea- 
surement. The initialization section prepares \+)\iljg{Bx, Bz)), and the probing sec- 
tion implements the global evolution U{t) approximately by decomposing it into 

^-iTB^{al+al)/2^-iTala'^^^-iTSa'^^al^-iTea°al^~iTB^{al+al)^-iTB^{al+al^ In the eXpCrimCUtS, WC 

measured the overlap L for a transverse field strength of = 0.1, coupling strengths of 
e = 0.2 and 0.3, and a range of longitudinal fields, —2 < B^ < 2. The evolution time was 
set to r = 1.6. The approximations used in the implementation of U{t) reduce the fidelity 
by less than 1.4%. 

The experimental results are shown as Figure [S, where the experimentally measured 
overlaps L are marked by "*" and "x" for e = 0.2 and 0.3, respectively. The experimental 
data are fitted to qLq, where Lq denotes the corresponding theoretical result. The best 
agreement was obtained for a = 0.84 and 0.77, respectively; the corresponding functions are 
shown as the dark and light curves. Obviously the critical points are correctly identified 
by the minima of L, indicating increased sensitivity of the ground state to the perturbation 
by the probe qubit. The differences between the theoretical and experimental values are 
mainly caused by imperfections of the radio frequency pulses, inhomogeneities of magnetic 
fields and decoherence. 

Discussion and Conclusion. — In conclusion, we have shown that a probe qubit can be 
used to detect quantum critical points. It is first placed into a superposition state and then 
coupled to the system undergoing the QPT. When the two eigenstates become correlated 
to two different phases, the superposition decoheres. The loss of coherence is thus a direct 
measure of the QPT. 

We have applied this procedure to two types of QPTs, choosing the couplings between 
the probe and the system in such a way that the two states of the probe induce slightly 
different values of the control parameter [12I, 17|. No details have to be known about the 



phases on the two sides of the phase transition. Only one qubit is measured for the detection 
of the critical points, independent of the size of the simulated quantum system. Hence this 
method scales very favorably with the size of the system [isl, . Theoretical results indicate 
that the overlap L remains a useful measure for larger systems in Ising and XY spin chains 
[1^ . For the more complex quantum phase transitions where many states are close to the 
ground state (e. g. spin glass), our fidelity method seems to work, although the details are 



still being worked out 19 1 
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In the present example, the probe qubit was coupled to all system qubits in a symmetric 
way. For other systems, the type of coupling required may depend on the system Hamiltonian 
and the nature of the phases on both sides of the QPT. While a full discussion of this issue 
is far beyond the scope of this letter, we expect that if the phase change involves delocalized 
states (e.g. spin waves), a single coupling between the probe qubit and one of the system 
qubits should be sufficient to detect the phase transition [2^. On the other hand, if the 
changes at the QPTs are local, a larger number of couplings or probe qubits may be required. 
In the extreme case, where critical points separate purely local changes, it may be necessary 
to couple the probe qubit to every system qubit or to implement couplings from a single 
probe qubit to all system qubits. Even in this worst case scenario, the number of probe qubits 
(or operations) only scales linearly with the size of the system; this should be contrasted 
to the readout by quantum state tomography, where the number of measurements increases 
exponentially with the system size. In future work, we plan to apply this type of analysis to 
the study of different types of phase transition, including quantum chaos 2l|. Furthermore, 
it should be possible to use this approach for the characterization of decoherence and 
errors that occur during quantum information processing 22|. 

We thank Prof. J.-F. Du for helpful discussions. This work is supported by the Alexander 
von Humboldt Foundation, the DFG through Su 192/19-1, and the GraduiertenkoUeg No. 
726. 
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FIG. 1: (Color online) (a) The energy levels of the system, (b) Concurrence (thick line) and overlap 
L (thin) of the ground state. 
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FIG. 2: (a) Quantum network for measuring L. W denotes the Walsh-Hadamard transform, and 
(7+ = {(Tx + i(Ty)/2. The controlled operations Pq and Pi denote the evolutions for preparing 
\'4'g{Bzfi)) and \iljg{Bz^i)) , if qubit is in state |0) or |1), respectively, (b) Chemical structure of 
Diethyl fluoromalonate. The three qubits are marked by the dashed oval, (c) Pulse sequence for 
measuring the state overlap. The narrow unfilled rectangles denote -7r/2 pulses, and the wide ones 
denote vr pulses. The striped rectangles denote 7r/4 pulses. The directions along which the pulses 
are applied are denoted by =ba; and zty. The durations of the pulses are so short that they can be 
ignored. 
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Magnetic field B. 



FIG. 3: Theoretical (line) and measured (asterisks) overlap L for the level-crossing case. Three 
NMR spectra illustrate the signals corresponding to Bz = —1.5, 0, and +1, respectively. 
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FIG. 4: (Color online) Quantum circuit (a) and pulse sequence (b) for measuring the overlap L 
for the avoided level-crossing case. P denotes the preparation of \+)\^g{Bx, B^)) from |000). In 
the initialization section, the width of the filled pulse applied to HI is — |, and the width of the 
filled pulses applied to F2, from left to right, are (vr — a)/2, (a — /3)/2, and (vr — /?)/2, respectively, 
with tan(a/2) = — -v/2ci/c+, tan(/?/2) = — c+/\/2co, tan(6'/2) = sin(/?/2)/ cos(a/2). cq, c+ and ci 
denote the amplitudes of |00), |</>+) and in \'ipg{Bx, Bz)) ■ In the probing section, the rectangles 
filled by heavy and light color denote the pulses with width tB^ and 2tBz, respectively, and 
= + 7^)' d2 = ^(^ + and ds = f + p^^)- 
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FIG. 5: (Color online) Experimental overlap L for e = 0.2 marked by "*" and e = 0.3 marked by 
"x". The experimental data are fitted to oLq, and yielded a = 0.84 and 0.77, respectively, shown 
as the dark and light curves, where Lq denotes the corresponding theoretical result. 
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